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PREFACE 


The  work  described  in  this  report  was  authorized  under  Project  No.  206023,  Non- 
Medical  CB  Defense.  This  work  was  started  in  October  2000  and  was  completed  in  February 


2003. 


Calculations  were  completed  using  the  Computing  Facilities  at  the  DoD  Major 
Shared  Resource  Center,  U.S.  Army  Research  Laboratory,  Aberdeen  Proving  Ground,  MD. 

The  use  of  either  trade  or  manufacturers’  names  in  this  report  does  not  constitute 
an  official  endorsement  of  any  commercial  products.  This  report  may  not  be  cited  for  purposes 
of  advertisement. 

This  report  has  been  approved  for  public  release.  Registered  users  should  request 
additional  copies  from  the  Defense  Technical  Information  Center;  unregistered  users  should 
direct  such  requests  to  the  National  Technical  Informational  Service. 
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HYDROLYSIS  OF  PHOSPHORUS  ESTERS: 
A  COMPUTATIONAL  STUDY 


1.  INTRODUCTION 

Computational  chemistry  is  a  valuable  tool  for  (1)  increasing  our  understanding  of 
chemical  phenomena,  (2)  elucidating  the  fundamental  principals  underlying  the  behavior  of  a 
diverse  set  of  chemical  compounds  under  varying  conditions,  and  (3)  reducing  the  number  of 
potential  chemical  candidates  for  a  given  study  to  a  manageable  set  of  characteristic  compounds 
that  can  be  examined  in  greater  detail. 

Organophosphorus  (OP)  nerve  agents  are  reactive  molecules  that  respond  to 
nucleophilic  attack.  Their  toxicity  results  from  inhibition  of  acetylcholinesterase  (AChE)  and  the 
subsequent  uncontrolled  accumulation  of  acetylcholine  at  the  neural  synapses.  The  germinal 
event  is  the  nucleophilic  attack  on  the  phosphorus  center  by  the  hydroxyl  moiety  in  a  serine 
residue  at  the  active  site  of  the  enzyme.  A  covalent  bond  is  formed  between  the  enzyme  and  the 
agent  with  the  simultaneous  dissociation  of  the  leaving  group  (i.e.,  F",  RS\  NC‘,  etc.). 

The  effectiveness  of  various  compounds  as  AChE  inhibitors  depends  upon  steric 

interactions  and  reactivity.  Compounds  that  are  too  reactive  frequently  hydrolyze  before 

reaching  the  critical  target  and  thereby  become  inactive.  Compounds  that  are  too  stable  might 

reach  the  target  and  bind  reversibly  to  AChE  but  be  resistant  to  nucleophilic  attack.  Quantum 

calculations  at  the  semiempirical  level  in  conjunction  with  mouse  lethality  data  indicated  that  the 

phosphonates  (i.e.,  G  agents)  have  optimal  reactivity  for  maximal  toxicity;  whereas,  the 
•  1,2 
phosphinates  are  too  reactive  and  phosphates  are  too  stable.  Also,  phosphonofluoridates  have 

optimal  reactivity;  whereas,  the  phosphonochloridates  are  too  reactive  for  CW  agents  but  are 

extremely  useful  as  reactive  intermediates  during  production.  The  optimal  reactivity  for  the 

leaving  group  depends  on  the  electronic  structure  of  the  entire  molecule.  In  the  mustard  series, 

the  chlorides  appear  to  be  optimal,  but  the  fluorides  are  too  stable  to  allow  formation  of  the 

cyclic  sulfonium  ions  at  effective  rates  under  ambient  conditions. 

The  reactivity  of  agents  is  also  a  critical  contributor  to  the  effectiveness  of 
decontamination  and  other  defensive  measures.  The  developers  must  select  decontaminants  that 

3 

react  within  a  few  minutes  to  form  nontoxic  materials.  Also,  it  is  required  that  the  chemicals 
not  destroy  the  equipment,  damage  the  environment,  or  injure  the  soldier.  Therefore,  selecting 
the  decontaminant  with  the  optimal  reactivity  is  crucial.  Most  of  the  decontaminants  are 
oxidizing  agents  that  convert  the  agents  into  a  less  reactive  species  or  nucleophiles  that  remove 

4 

the  leaving  group  and  make  the  agent  less  reactive. 

The  stability  of  OP  compounds  is  dependent  on  the  leaving  group  and  the 
substituents  attached  to  the  phosphorus  atom.  A  poor  leaving  group  can  either  slow  the  reaction 
rate  or  render  the  OP  compound  unreactive.  Also,  the  electronic  characteristics  of  substituents 
attached  to  the  phosphorus  atom  play  an  important  role  in  stabilizing  or  destabilizing  the 
transition  state  (TS)  and  therefore  affect  the  reaction  rates.  This  is  observed  in  the  relative 
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reactivities  of  phosphinates  (two  P-alkyl  bonds),  phosphonates  (one  P-alkyl  bond  and  on  P- 

5-7 

alkoxy  bond),  and  phosphates  (two  P-alkoxy  bonds).  In  other  words,  the  activation  energy  for 
the  phosphinate  will  be  lower  than  the  phosphonate,  which  will  be  lower  than  the  phosphate. 
Because  the  P-alkyl  bond  increases  reactivity  (relative  to  a  P-alkoxy  bond),  phosphinates  are  the 
most  and  phosphate  the  least  reactive  of  the  three  classes  of  phosphorus  esters. 


2.  RESEARCH  PLAN 

To  better  understand  the  mechanisms  for  the  hydrolysis  of  CW  agents,  three  types 
of  OP  compounds;  phosphinates,  phosphonates,  and  phosphates  have  been  investigated  using  an 
anionic  (hydroxide)  and  neutral  (water)  reaction  pathway  schemes.  In  the  anionic  pathway,  the 
OP  compounds  have  been  investigated  in  the  gas  phase  and  in  solution  (using  the  SCI-PCM 

solvation  method).  In  the  neutral  reaction  pathway,  which  is  similar  to  our  previous  report  on  the 

8 

methanolysis  of  the  phosphinate,  gas  phase  results  augmented  by  explicit  solvation  models  are 
presented  instead  of  SCI-PCM  calculations  since  the  latter  exhibited  significant  variational 
instability.  This  study  is  not  only  intended  to  address  the  hydrolysis  mechanism,  but  also  to  give 
insight  relevant  to  investigations  of  the  phosphorylation  mechanisms  in  the  active  site  of  AChE. 
As  an  example,  in  the  active  site  of  AChE,  the  catalytic  efficacy  is  achieved  through 
simultaneous  participation  of  several  residues  including  the  catalytic  triad  and  oxyanion  hole. 

The  oxyanion  hole  stabilizes  the  negative  charge  on  the  phosphoryl  bond  and  in  some  ways 
mimics  how  water  molecules  solvate  the  OP  compound  during  hydrolysis. 


3.  COMPUTATIONAL  PROCEDURES 

All  ab  initio  molecular  orbital  calculations  were  performed  using  Gaussian 

9 

Incorporated’s  G94  (anionic  mechanism)  and  G98  (for  all  other  calculations)  programs.  Initial 

geometries  for  all  structures  were  obtained  from  low-level  calculations  at  the  RHF/3-21G  level 

of  theory.  All  structures  for  the  anionic  reaction  pathways  were  then  fully  optimized  using  the 

B3LYP  exchange  and  correlation  functionals1011  with  the  6-31+G(2d)  basis  set.  A  method- 

dependence  study  was  carried  out  on  the  neutral  reaction  methodology  (4-membered  transition 

state),  by  investigating  the  complete  reaction  pathway  for  the  phosphinate  system  using  the  RHF 

level  of  theory  with  the  6-3 1  lG(2d,2p)  basis  set,  the  B3LYP  exchange  and  correlation 
10,11 

functionals  with  the  6-31 1+G(2d,2p)  basis  set,  and  the  second-order  Moller-Plesset 

12,13 

perturbation  theory  (MP2)  with  the  6-3 1 1+G(2d,2p)  basis  set.  After  analyzing  the  data,  it 

was  observed  that  the  DFT  gave  comparable  structural  and  energetic  data  relative  to  the  more 

computationally  demanding  MP2  level  of  theory,  therefore  the  pathways  for  the  corresponding 

phosphonate  and  phosphate  were  only  investigated  using  the  B3LYP/6-31 1+G(2d,2p)  level  of 

theory.  In  the  neutral  reaction  methodology  (6-membered  transition  state),  the  complete  reaction 

pathway  for  all  three  OP  compounds  was  investigated  using  the  B3LYP  exchange  and 

10,11 

correlation  functionals  with  the  6-3 1 1+G(2d,2p)  basis  set.  For  both  neutral  reaction 
pathways,  the  SCI-PCM  SCRF  solvation  model  could  not  be  employed  successfully  due  to 
convergence  difficulties  associated  with  the  proton  transfer  in  these  systems.  All  minimizations 

14,15 

were  carried  out  using  the  Bemy  algorithm  and  the  default  parameters  were  used  for  the 
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integral  cut-off  and  minimization  convergence  criteria.  Vibrational  frequencies  were  calculated 
for  all  critical  points  to  (1)  obtain  the  zero  point  and  Gibbs  free  energies  and  (2)  insure  that  the 
transition  states  had  only  one  imaginary  frequency  while  all  minima  had  zero  imaginary 

16,17 

frequencies.  Intrinsic  Reaction  Coordinate  (IRC)  calculations  were  performed  on  the  gas 
phase  transition  states  for  the  anionic  pathways  and  the  phosphinate  neutral  reaction  pathway  to 
confirm  that  the  transition  states  did  in  fact  connect  their  corresponding  minima. 

For  the  inclusion  of  solvation  on  the  anionic  systems,  the  Onsager  SCRF 

18,19 

model  was  initially  attempted,  with  results  no  significant  differences  from  the  isolated 

20,21  22 

results.  Due  to  an  unavailability  of  second  derivatives  within  the  PCM  and  IPCM 
methods,  the  self-consistent  isodensity  polarized  continuum  method  (SCI-PCM)22  implemented 
in  GAUSSIAN-94  was  used  instead.  The  base  level  of  theory  was  the  same  as  the  gas  phase 
calculations.  A  dielectric  constant  of  78.0  (approximately  that  of  room  temperature  H2O)  was 
used  for  the  solvent.  A  value  for  the  isodensity  surface  of  0.0002  au  was  used  for  all  systems, 
thus  eliminating  the  convergence  criterion  failures  observed  with  the  recommended  0.0004  au 
value  on  the  larger  systems.  The  internal  special  grid  of  974  was  used  to  define  the  total  number 
of  phi  and  Gauss-Legendre  theta  values  for  numerical  surface  integration. 


4.  RESULTS/DISCUSSION 

4.1  Anionic  Reaction  Pathway. 

A  general  reaction  schematic  for  all  three  OP  anionic  systems  with  the  significant 
atoms  numbered  is  presented  in  Scheme  1.  The  hydroxide  anion  approaches  the  phosphorus 
atom  anti  to  the  fluoride-leaving  group  forming  a  trigonal  bipyramidal  (TB)  intermediate  with 

Scheme  1 
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the  nucleophile  and  leaving  group  in  the  axial  positions  before  forming  the  products.  For 
simplicity  in  discussing  results  on  the  anionic  system,  the  gas  phase  and  SCI-PCM  results  for 
each  of  the  three  anionic  OP  systems  are  presented  together  in  the  text  and  corresponding 
figures.  In  the  reaction  pathway  figures,  both  potential  energy  surfaces  (PES)  are  displayed,  with 
the  gas  phase  optimized  structures  depicted  relative  to  their  positions  along  the  reaction 
coordinate.  The  reactants  and  products  at  infinite  separation  are  not  depicted. 
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The  isolated  and  solvated  anionic  PESs  for  the  hydrolysis  of  the  phosphinate, 
dimethylphosphinic  fluoride,  are  depicted  in  Figure  1  with  the  corresponding  energetic  and 
structural  results  in  Table  1.  In  the  gas  phase  free  energy  profile  in  Figure  1  (dashed  line),  the 
transition  states  for  the  incoming  nucleophile  and  the  departing  leaving  group  are  below  the 
reactants  and  products  at  infinite  separation.  This  trend  is  expected  because  the  anions  are 
very  unstable  in  the  absence  of  solvent  stabilization  effects.  Therefore,  a  correlation  between  the 
reaction  rates  of  the  three  OP  systems  can  not  be  deduced  from  the  gas  phase  results,  however, 
the  gas  phase  structures  are  used  to  help  understand  the  complete  reaction  pathway  and  are 
helpful  in  yielding  starting  geometries  for  the  solvation  methods. 


Figure  1.  Free  Energies  for  the  Anionic  Paths  for  the  Hydrolysis 
of  Dimethylphosphinic  Fluoride 
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Free  energy  profiles  for  the  gas  phase  and  SCI-PCM  anionic  reactions.  Gas  phase  optimized  geometries  are 
shown  relative  to  their  position  along  the  reaction  coordinate.  Methyl  hydrogens  omitted  for  clarity. 
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Table  1.  Hydroxide  with  Dimethylphosphinic  Fluoride  —  Structural  and  Energetic  Results* 
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phase  and  in  solution  (SCI-PCM)  respectively.  The  relative  free  energies  (in  kcal/mol)  are  relative  to  the  lowest 
energy  structure.  Bond  lengths  are  in  angstroms  and  angles  in  degrees. 


The  hydroxide  anion  forms  a  stable  hydrogen  bonded  (HB)  complex  2  with  the 
methyl  protons  of  the  phosphinate  in  an  orientation  opposite  of  the  fluoride-leaving  group.  Even 
though  the  hydroxide  oxygen,  03  (Scheme  1)  is  3. 5 A  from  the  phosphorus  atom,  long-range 
charge  stabilization  occurs  through  the  methyl  groups  via  the  slight  lengthening  of  the  PI -03  and 
P1-F4  bonds  in  both  systems  (Table  1).  The  geometries  for  the  HB  complex  2  in  both  systems 
are  very  similar  with  the  exception  of  a  .10°  angle  difference  for  the  incoming  nucleophile 
(Table  1:  rc03-Pl-02  and  7i03-Pl-F4).  The  hydroxide  oxygen  then  forms  a  covalent  bond  with 
the  phosphorus  atom  to  form  the  stable  TB  structure,  4,  via  TS  3.  The  geometries  for  TS  3  in 
both  systems  are  very  similar  with  a  slightly  longer  PI -03  partial  bond  distance  arising  from  the 
inclusion  of  solvation  effects.  In  TS  3,  the  negative  charge  is  distributed  to  the  P1-F4  bond  and 
not  the  Pl-02  bond  (Table  1).  Looking  at  the  SCI-PCM  PES  in  Figure  1,  an  energy  barrier  of 
4.7  kcal/mol  from  reactants  to  structure  4  is  now  observed  via  TS  3.  The  TB  structure  4  is 
slightly  distorted  with  an  angle  of  167°  in  the  axial  substituents  (7i03-Pl-F4  in  Table  1).  A 
difference  in  the  Pl-03  and  P1-F4  bond  lengths  is  also  observed  between  the  two  methods  with 
the  solvated  calculations  being  better  able  to  stabilize  the  charge  on  the  leaving  group.  Also  in 
structure  4,  the  single  bond  character  associated  with  the  negative  charge  located  on  the 
phosphinyl  oxygen  is  observed.  The  fluoride  ion  departs  from  structure  4  via  TS  5  with  a  7.0 
and  1.6  kcal/mol  energy  barrier,  respectively  (Table  1).  Here  it  is  observed  that  the  existence  of 
the  stable  pentacoordinate  species  in  the  gas  phase  is  minimized  in  solution.  HB  complex  6  was 
only  treated  in  the  gas  phase  due  to  convergence  criterion  problems  arising  from  the  frequency 
calculation  of  the  SCI-PCM  method.  An  overall  exothermic  reaction  pathway  is  observed  for  the 
SCI-PCM  PES  with  an  initial  energy  barrier  of  4.7  kcal/mol. 

The  isolated  and  solvated  anionic  PESs  for  the  hydrolysis  of  the  phosphonate, 
methyl  methylphosphonofluoridate,  are  depicted  in  Figure  2  with  the  corresponding  energetic 
and  structural  results  reported  in  Table  2.  The  reaction  pathways  are  very  similar  to  the 
phosphinate  with  the  exception  of  an  additional  energy  barrier,  TS  5,  for  the  rotation  of  the 
methoxy  substituent.  The  hydroxide  forms  a  similar  HB  complex  2  to  the  protons  on  the  methyl 
groups  as  in  the  phosphinate  system.  The  lengthening  of  the  P1-F4  bond  is  not  as  significant  as 
in  the  phosphinate  system,  which  is  due  to  a  slight  lengthening  of  the  C-O  methoxy  bond  (not 
shown  in  Table  2).  The  stable  TB  structure  4  is  formed  via  TS  3.  In  the  SCI-PCM  method  an 
energy  barrier  of  6.21  kcal/mol  is  observed  from  reactants  to  structure  4.  Again  there  is  a 
distortion  of  a  true  TB  geometry  with  a  7t03-Pl-F4  of  168.8  and  169.7E,  respectively.  To 
prepare  for  the  departure  of  the  fluoride-leaving  group,  the  methoxy  group  needs  to  be  on  the 
same  side  as  the  leaving  group.  This  is  accomplished  through  a  small  methoxy  rotation  of  2.76 
and  2.74  kcal/mol,  respectively,  via  TS  5  to  form  the  other  TB  structure  6.  Structure  6  is  more 
stable  than  structure  4  by  1.01  and  0.57  kcal/mol,  respectively.  Therefore,  there  is  a  preference 
for  this  rotation  to  occur.  With  the  methoxy  on  the  same  side  as  the  leaving  group,  the  fluoride 
departs  through  a  5.54  and  3.21  kcal/mol  energy  barrier  via  TS  7,  respectively.  Again,  due  to 
convergence  criterion  problems  arising  from  the  frequency  calculation  of  the  SCI-PCM  method, 
HB  complex  8  was  only  found  in  the  gas  phase.  The  isolated  and  solvated  anionic  PESs  for  the 
hydrolysis  of  the  phosphate,  dimethyl  phosphorofluoridate,  in  the  gas  phase  and  in  solution  are 
depicted  in  Figure  3  with  the  corresponding  energetic  and  structural  results  in  Table  3.  Again, 
the  reaction  pathway  is  similar  to  the  two  previous  pathways  with  the  exception  of  one  more 
methoxy  rotation.  In  HB  complex  2,  there  is  no  lengthening  of  the  Pl-02  or  P1-F4  bonds  in 
either  system.  In  fact,  in  the  SCI-PCM  method  there  is  a  slight  shortening  of  the  P1-F4  bond 
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(Table  3).  In  this  system,  the  distribution  of  the  negative  charge  is  placed  equally  on  the  C-0 
bonds  of  the  methoxy  groups  (not  shown  in  Table  3).  In  the  SCI-PCM  method,  there  is  a 
6.66  kcal/mol  energy  barrier  from  reactants  to  the  first  TB  structure  4  via  TS  3.  Again,  the 
geometries  for  TS  3  are  similar  in  both  methods  with  the  exception  of  the  PI -03  distance,  which 
is  slightly  longer  with  the  inclusion  of  solvation.  As  depicted  in  Figure  3,  structure  4  has  both 
methoxy  groups  on  the  side  of  the  incoming  nucleophile.  To  aid  in  the  departure  of  the  leaving 
group,  both  methoxy  groups  need  to  be  on  the  same  side  as  the  fluoride-leaving  group.  This  is 
accomplished  (for  the  purpose  of  plotting  the  complete  reaction  pathway)  through  a  step-wise 
process  of  one  methoxy  rotation  at  a  time.  The  first  rotation,  as  depicted  in  Figure  3,  is  the  back 
methoxy  group  with  an  energy  barrier  of  2.73  and  3.11  kcal/mol,  respectively,  forming  TB 
structure  6.  The  second  rotation  is  the  front  methoxy  group  with  an  energy  barrier  of  2.74  and 
2.77  kcal/mol,  respectively,  forming  TB  structure  8.  As  with  the  phosphonate  system,  there  is  a 
thermodynamic  preference  for  these  rotations  to  occur.  Structure  8  is  more  stable  than  structure 
4  by  2.21  and  0.87  kcal/mol,  respectively  (Table  3).  From  structure  8,  the  fluoride  departs  with  a 
7.86  and  5.93  kcal/mol  energy  barrier  via  TS  9.  HB  structure  10  could  not  be  located  with  the 
SCI-PCM  solvation  method.  It  is  also  worth  mentioning  that  throughout  the  reaction  pathway, 
even  through  the  methoxy  rotations,  there  is  a  trend  of  the  PI -03  shortening  and  the  P1-F4 
lengthening  (Table  3). 


Figure  2.  Free  Energies  for  the  Anionic  Paths  for  the 

Hydrolysis  of  Methyl  Methylphosphonofluoridate 


Gas  phase  optimized  geometries  are  shown  relative  to  their  position  along  the  reaction  coordinate.  Methyl 
hydrogens  omitted  for  clarity. 
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Table  2.  Hydroxide  with  Methyl  Methylphosphonofluoridate  —  Structural  and  Energetic  Results* 
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*The  two  values  shown  for  each  parameter  are  optimized  using  the  B3LYP/6-31+G(2d)  level  of  theory  in  the  gas 
phase  and  in  solution  (SCI-PCM),  respectively.  The  relative  free  energies  (in  keal/mol)  are  relative  to  the  lowest 
energy  structure.  Bond  lengths  are  in  angstroms  and  angles  in  degrees. 


Table  3.  Hydroxide  with  Dimethyl  Phosphorofluoridate  —  Structural  and  Energetic  Results* 
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*The  two  values  shown  for  each  parameter  are  optimized  using  the  B3LYP/6-31+G(2d)  level  of  theory 
in  the  gas  phase  and  in  solution  (SCI-PCM),  respectively.  The  relative  free  energies  (in  kcal/mol)  are 
relative  to  the  lowest  energy  structure.  Bond  lengths  are  in  angstroms  and  angles  in  degrees. 


Figure  3.  Free  Energies  for  the  Anionic  Paths  for  the 
Hydrolysis  of  Dimethyl  Phosphorofluoridate 


Free  energy  profiles  for  the  gas  phase  and  SCI-PCM  anionic  reaction  pathways  for  the  hydrolysis  of  dimethyl 
phosphorofluoridate.  Gas  phase  optimized  geometries  are  shown  relative  to  their  position  along  the  reaction 
coordinate.  Methyl  hydrogens  omitted  for  clarity. 


As  can  be  observed,  the  trends  of  phosphinates  hydrolyzing  faster  than 
phosphonates  which  hydrolyze  faster  than  phosphates  are  reproduced  with  the  SCI-PCM  method 
with  energy  barriers  of  4.67,  6.03,  and  6.66  kcal/mol,  respectively.  However,  the  energy  barriers 
are  significantly  too  low  and  close  together.  These  barriers  are  a  couple  of  kcal/mol  below  the 

5,6 

respective  chlorides  ,  which  are  more  reactive  than  the  fluorides.  Therefore,  the  SCI-PCM 
method  predicts  the  correct  order  (i.e.,  qualitatively  correct)  but  inadequately  describes  the  rates 
of  reaction  (i.e.,  quantitatively  incorrect).  The  SCI-PCM  method  does,  however,  stabilize  the 
negative  charge  that  is  observed  in  the  differences  of  the  PI -03,  PI -02  and  P1-F4  bond  lengths, 
and  it  provides  a  superior  representation  of  the  overall  reaction  pathway  relative  to  the  gas  phase 
results.  Apparently  a  more  rigorous  solvation  method  needs  to  be  implemented  to  accurately 
describe  the  energy  barriers,  however. 
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4.2 


Neutral  Reaction  Pathway  -  4-Membered  Transition  State. 


In  an  attempt  to  systematically  investigate  the  hydrolysis  of  the  three  OP 
compounds  using  discrete  water  molecules  to  mimic  solvation  effects  arising  from  a  shell  of 
water,  a  methodology  starting  with  only  one  water  molecule  was  undertaken.  Such  methodology 
using  a  single  water  molecule  leads  to  a  somewhat  strained  4-membered  TS  that  is  closely 

analogous  to  a  reaction  that  has  been  previously  demonstrated.  Whereby,  the  methanolysis  of 

8 

dimethylphosphinic  fluoride  mimics  the  serine  group  of  AChE.  A  schematic  for  the  neutral 
reaction  pathway  is  depicted  in  Scheme  2  with  the  significant  atoms  numbered.  The  nucleophile 
03  oxygen  approaches  the  phosphorus  atom  anti  to  the  fluorine-leaving  group  as  in  the  anionic 
reaction  pathway.  However  in  this  methodology,  the  H5  proton  from  the  nucleophile  transfers  to 
the  phosphonyl  02  oxygen  (stabilizing  the  negative  phosphoryl  group)  simultaneous  with  the 
03-PI  formation.  The  H5  proton  then  departs  with  the  fluoride-leaving  group,  forming  the 
phosphate  and  hydrogen  fluoride  molecule.  The  basis  for  using  this  method  is  in  the  use  of  the 
H5  proton  from  the  water  molecule  to  stabilize  the  negative  charge  in  place  of  the  effective 
solvation  shell  model  arising  from  the  SCRF  method.  This  explicit  molecule  stabilization  can 
also  be  compared  to  the  role  of  the  oxyanion  hole  in  the  active  site  of  AChE  during  the 
phosphorylation  of  the  serine  hydroxyl  group.  If  the  role  of  solvation  methods  is  to 
reduce/stabilize  charges,  then  this  neutral  reaction  methodology  represents  an  “ultimate” 
solvation  model.  It  needs  to  be  stressed,  however,  that  the  above  neutral  reaction  pathway  is  a 
general  methodology  and  should  be  viewed  as  a  first  approximation  in  a  series  of  models 
including  increasingly  large  numbers  of  water  molecules,  and  is  intended  to  obtain  a 
quantitatively  reasonable  representation  of  a  full  solvation  sphere  at  a  reasonable  computational 
cost. 


Scheme  2 


IT 


H5 


'03 


+ 


02 

II 

R'yP  1' 


R 


"F4 


H5. 


H 


02 


03— PI — F4 
R  R' 


02 

II 

Plu„ 


H03' 


R 


+ 


T 


Phosphinate: 

Phosphonate: 

Phosphate: 


R=R'=Me 
R=Me,  R-OMe 
R=R-OMe 


The  structural  and  energetic  results  for  the  complete  gas  phase  one-water  neutral 
reaction  pathways  for  the  hydrolysis  of  dimethylphosphinic  fluoride  using  the  RHF/6- 
3 1 1  G(2d,2p),  B3LYP/6-3 1 1+G(2d,2p),  and  MP2/6-31 1+G(2d,2p)  levels  of  theory  are 
consolidated  in  Table  4.  For  simplicity  in  discussing  results  in  the  text  of  the  neutral  system,  the 
first,  second  and  third  sets  of  results  reported  are  at  the  RHF,  B3LYP,  and  MP2  levels  of  theory, 
respectively.  In  Figure  4,  the  PESs  of  the  RHF  and  B3LYP  levels  of  theory  are  depicted  together 
with  the  B3LYP  gas-phase  optimized  structures  relative  to  their  positions  along  the  reaction 
coordinate.  The  MP2/6-3 1 1+G(2d,2p)  energetics  are  not  depicted  in  Figure  4  due  to  very  close 
overlap  with  the  DFT  level.  The  reactants  and  products  at  infinite  separation  are  not  depicted. 
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Table  4.  One-Water  Hydrolysis  of  Dimethylphosphinic  Fluoride  —  Structural  and  Energetic  Results* 


w 


tO 

VO 

IT) 

CN 

•O 

CN 

T— H 

NO 

VO 

l 

o 

VO 

o 

oo 

t— H 

t— H 

VO 

tO 

tO 

ON 

O 

ON 

On 

O 

1 

CO 

ON 

l> 

ON 

VO 

00 

VO 

00 

o 

1—1 

r- 

r- 

00 

td 

Tl* 

**fr 

CN 

ON 

ON 

o 

r-H 

ri 

o 

r- 

t > 

r- 

r- 

r-H 

r- 

^H 

r- 

r- 

t— H 

r— H 

T— H 

L- 

r-H 

VO 

VO 

i— 1 

i—l 

VO 

T— H 

t— H 

VO 

omvu'i'tw 
Tt  </■>  O  «/")  Tj-  O 

ri  ri  rt  d  d  ri 


m 

r> 

Tf 

VO 

CN 

CN 

00 

o 

tO 

<0 

00 

tO 

CN 

ON 

^H 

vo 

tO 

ON 

o 

ON 

00 

00 

On 

<N 

o 

tO 

cq 

vo 

o 

o 

r— H 

vo 

ON 

cd 

00 

vd 

vd 

vd 

00 

00 

On 

vd 

to 

to 

od 

ON 

o 

o 

t-H 

00 

ON 

ON 

On 

CO 

CO 

00 

00 

00 

00 

CO 

00 

00 

vo 

VO 

r- 

to 

< 


m  rs  os  on 
vq  n  N 

f-H  00  — I  ■— 1 1  fT 

VO-lft  VO  h  >• 


>0 

ON 

V*) 

vo 

>0 

ON 

V) 

CC 

^H 

tO 

m 

00 

vq 

p 

00 

t-h 

rc 

vq 

rn 

t-H 

CN 

cn 

cd 

rd 

ON 

ON 

od 

r-H 

O* 

o 

(— s 

O 

o 

ON 

ON 

On 

ON 

On 

ON 

00 

00 

00 

t-H 

^H 

«; 


£ 


so  o  os  csi  h- 
p  (N  co  r|  p 
«-i  r  i  f  i  rn  nr] 


(Si 

o 

1—1 

o 

to 

ON 

o 

00 

to 

m 

vo 

o 

00 

vo 

CO 

(N 

Tt* 

o 

t^ 

to 

ON 

00 

to 

o 

ON 

tO 

ON 

00 

vo 

cn 

r-H 

On 

CO 

ro 

t-H 

o 

vo 

<N 

r—H 

CO 

1- * 

00 

t-H 

«o 

vo 

tO 

»o 

to 

to 

to 

vq 

vo 

vq 

vq 

vq 

p 

fO 

(N 

fTi 

Tt; 

rn 

r—H 

T-l 

i—H 

t— i 

1— 1 

t-H 

i-H 

1-1 

t-H 

i— H 

r-H 

i— H 

i— ^ 

1—1 

1-1 

r-H 

oi 

oi 

oi 

rd 

m 

cd 

iz 


to 

ro 

i— i 

(N 

vo 

o 

to 

TT 

r-H 

vo 

to 

00 

vo 

<N 

i-H 

oo 

to 

o 

t" 

CN 

t-H 

o 

< 

o 

i-H 

o 

r- 

On 

m 

r- 

t-H 

o 

to 

00 

ON 

ON 

CNI 

00 

<N 

t-H 

On 

m 

(Nl 

i 

rn 

rq 

<N 

CN 

i-H 

vq 

p 

p 

vq 

vq 

vq 

to 

vq 

vo 

to 

vq 

vq 

to 

vq 

vq 

r-H 

IP 

rd 

rd 

rd 

(N 

cd 

cd 

r-M 

T-H 

i-H 

1— H 

t-H 

T-H 

i—H 

t-H 

r-H 

1— H 

t-H 

(Nl 

ON 

o 

-»H- 

00 

CN 

t-H 

to 

•o 

to 

r- 

m 

ON 

ON 

o 

00 

T-H 

On 

(O 

Tf 

o 

o 

t— H 

to 

o 

(Tl 

F* 

00 

<N 

to 

to 

On 

<N 

CN 

o 

m 

ON 

m 

to 

m 

r- 

vo 

vo 

ON 

ON 

to 

00 

00 

1 

Tf 

to 

to 

to 

p 

vo 

vq 

p 

vo 

to 

p 

vq 

to 

to 

to 

rh 

r-H 

P 

i— H 

r— H 

t-H 

t-H 

r-H 

^H 

t-H 

t-H 

r-H 

T— H 

t-H 

rH 

t-H 

r-H 

1—1 

rH 

t-H 

rH 

0  0  0  0*00 
O  O  O  IT)  o  r- 


— 1  H- 
rf  IT) 


00 

On 

to 

r-H 

vo 

ro 

CN 

00 

to 

to 

CN 

CN 

T^H 

CN 

ON 

On 

00 

00 

ro 

r-I 

00 

to 

^H 

Tj- 

_ _ j 

oo 

ON 

to 

r—H 

ON 

vd 

ON 

cd 

t-H 

CN 

CN 

CN 

CN 

i— H 

1 

CN 

CN 

CN 

,  co  n  n  h 
'  h  h  in  *h 


m  «-*  w> 

N  «  H 

vo  mi  co 


p 

o 


oq 

© 


CN 

t-H 

CN 

ON 

to 

o 

o 

o 

vo 

C- 

<o 

O 

to 

CN 

p 

o 

00 

to 

T—l 

t> 

to 

co 

o 

rH 

ON 

O 

r-H 

K 

vd 

vd 

cd 

cd 

vd 

cd 

CN 

ON 

CN 

r^j 

CN 

t-H 

CN 

CN 

CN 

CN 

t-H 

Os  o 
p  p 
o  © 


O  — I  CO  o  o  o 
0^000*0 
oo  K  os  ©  ©  © 


t-~  ■— <  Tt 
CO  OS  o 
©  © 


o 

rO 

00 

VO 

CN 

00 

VO 

00 

m 

to 

CN 

fTj 

t-H 

O 

rf 

00 

vo 

p 

p 

OO 

CN 

00 

rq 

rq 

ov 

00 

00 

00 

ON 

P 

t-H 

O 

G\ 

ON 

cd 

t-H 

T— H 

od 

to 

to 

^H 

r-H 

o 

o 

ON 

T— H 

<3 

o 

o 

CN 

CN 

Tt 

(O 

r^i 

T— H 

i— H 

1— H 

CN 

<N 

CN 

~4 

r-H 

T— H 

CN 

t-H 

CN 

r—H 

t-H 

Oi 


m 

00 

H 


vs 

oo 

H 


r 

oo 

H 


00 


Os 


18 


*The  values  for  each  parameter  are  optimized  using  the  RHF/6-31  lG(2d,2p),  B3LYP/6-311+G(2d,2p),  and  MP2/6- 
3 1 1+G(2d,2p)  levels  of  theory,  respectively.  The  relative  energies  (in  kcal/mol)  are  relative  to  the  lowest  energy  structure. 
E(z)  and  E(G)  are  the  zero-point  and  Gibbs  free  energies,  respectively.  Bond  lengths  are  in  angstroms  and  angles  in  degrees. 


kcal/mol 


Figure  4.  Relative  Free  Energies  for  the  One-Water  Neutral  Reaction  Pathway 
for  the  Hydrolysis  of  Dimethylphosphinic  Fluoride 
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Gas  phase  optimized  geometries  are  shown  relative  to  their  position  along  the  reaction  coordinate.  Methyl 
hydrogens  omitted  for  clarity. 


For  all  neutral  reaction  pathways  the  electronic,  zero-point  and  Gibbs  free 
energies  will  be  reported.  The  activation  energies  reported  in  the  text  will  be  the  Gibbs  tree 
energy  from  reactants  at  infinite  separation  to  TS  3  unless  otherwise  stated.  Where  appropriate, 
comparisons  of  the  electronic  and  zero-point  energies  from  the  hydrogen-bonded  structure  2  to 
TS  3  will  be  discussed. 

Structure  2  has  hydrogen  bonding  from  H5-02  and  from  03 -methyl  protons, 
which  aligns  the  water  molecule  anti  to  the  fluoride-leaving  group  (Figure  4).  In  TS  3,  the  H5 
proton  transfer  to  the  phosphinyl  oxygen  occurs  simultaneous  with  a  PI -03  bond  formation. 
Also  in  TS  3,  the  PI -02  bond  has  a  bond  order  of  1.5  with  a  1.52, 1.56,  and  1.56A  bond  lengths, 
respectively.  TB  structure  4  is  formed  via  TS  3  with  an  energy  barrier  of  44.4,  33.5,  and 
31.8-kcal/mol  energy,  respectively  (Table  4).  The  HF  level  of  theory  has  a  larger  activation 
energy  relative  to  the  other  two  levels  of  theory  as  expected,  however,  in  the  neutral  reaction 
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methodology,  the  TS  at  the  B3LYP  level  of  theory  is  approximately  2  kcal/mol  larger  than  the 

23 

MP2  level  of  theory,  which  is  in  contrast  with  other  proton  transfer  investigations.  When 
comparing  the  molecular  geometry  predictions  at  the  three  levels  of  theory,  the  B3LYP  and  MP2 
levels  lead  to  very  similar  minima  structures  with  only  minor  differences  in  the  TSs  (Table  4).  In 
structure  4,  the  7t03-Pl-F4  angle  of  178°  is  more  in  keeping  with  the  TB  geometry  than  was  the 
case  in  the  anionic  system.  For  the  proton  to  assist  in  the  fluoride  departure,  it  needs  to  be  on  the 
same  side  as  the  leaving  group.  This  is  accomplished  through  a  5.4,  4.8  and  3.6  kcal/mol 
rotation  around  the  phosphinyl  bond  to  form  TB  structure  6,  which  are  4.7,  4.8  and  4.1  kcal/mol 
more  stable  than  structure  4.  The  leaving  group  departs  via  TS  7  with  10.0-,  6.1-,  and 
7.0-kcal/mol  energy  barriers.  TS  7  has  simultaneous  proton  transfer  with  the  P1-F4  bond 
breakage  to  form  the  phosphate  and  hydrogen  fluoride,  structure  8.  Again,  it  is  important  to 
emphasize  that  the  proton  is  used  to  stabilize  the  charge  in  this  neutral  reaction.  In  a  real 
solvation  system,  several  water  molecules  would  stabilize  the  departing  fluoride  anion,  however, 
the  overall  mechanism  should  still  be  comparable. 

Due  to  the  similarities  of  the  structural  and  energetic  results  between  the  DFT  and 
MP2  levels  of  theory  of  the  phosphinate  (Table  4),  the  one- water  neutral  reaction  pathway  for 
the  hydrolysis  of  the  phosphonate  and  phosphate  were  carried  out  using  only  the  B3LYP/6- 
31 1+G(2d,2p)  level  of  theory  and  are  depicted  in  Figures  5  and  6,  respectively  with  the 
corresponding  energetic  and  structural  results  consolidated  in  Table  5.  The  B3LYP  data  for  the 
phosphinate  from  Table  4  is  also  included  in  Table  5  for  ease  of  comparison;  therefore,  the  three 
results  reported  in  the  text  will  correspond  to  the  phosphinate,  phosphonate  and  phosphate, 
respectively.  Due  to  the  number  of  possible  conformations  with  the  larger  OP  compounds,  the 
pathway  between  the  two  TB  structures  is  omitted.  The  reactants  and  products  at  infinite 
separation  are  not  depicted. 

It  deserves  mentioning  that  in  the  phosphonate  and  phosphate  systems,  the  lowest 
energy  conformation  found  is  depicted  and  discussed.  As  an  example;  for  the  first  TS,  the 
methoxy  group  oriented  on  the  same  side  of  the  fluorine  in  both  systems  (Figures  5  and  6)  is 
more  stable  than  the  methoxy  group  oriented  on  the  side  of  the  incoming  water(s)  (albeit  by  only 
a  very  small  amount). 

The  overall  PESs  for  the  phosphonate  and  phosphate  are  similar  to  the  PES  of  the 
phosphinate,  therefore;  only  noticeable  trends  will  be  contrasted.  In  comparing  the  PI -02,  Pl- 
03  and  P1-F4  bonds,  there  is  a  consistent  trend  of  the  phosphate  having  tighter/shorter  bonds 
than  the  phosphinate.  As  described  in  the  introduction,  this  correlates  well  to  the  impact  that 
substituents  have  on  the  relative  reactivities  of  the  three  classes  of  OP  compounds  studied.  It  is 
also  interesting  to  note  that  TB  structure  4  is  more  stable  than  TB  structure  6  only  for  the 
phosphate.  This  may  play  an  additional  role  to  slow  the  reactivity  of  the  phosphates.  In 
comparing  the  three  OP  systems  at  the  B3LYP/6-31 1+G(2d,2p)  level  of  theory,  the  free  energy 
barriers  are  33.5,  37.2  and  42.1  kcal/mol,  respectively. 
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Figure  5.  Relative  Free  Energies  for  the  One-Water  Neutral  Reaction  Pathway 
for  the  Hydrolysis  of  Methyl  Methylphosphonofluoridate 


methylphosphonofluoridate.  Gas  phase  optimized  geometries  are  shown  relative  to  their  position  along 
the  reaction  coordinate.  Methvl  hydrogens  omitted  for  claritv. 


Figure  6.  Relative  Free  Energies  for  the  One-Water  Neutral  Reaction  Pathway 
for  the  Hydrolysis  of  Dimethyl  Phosphorofluoridate 


phosphorofluoridate.  Gas  phase  optimized  geometries  are  shown  relative  to  their  position  along  the  reaction 
coordinate.  Methyl  hydrogens  omitted  for  clarity. 
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Table  5.  Structural  and  Energetic  Results  for  the  One- Water  Hydrolysis  of  Dimethylphosphinic  Fluoride,  Methyl 
Methylphosphonofluoridate,  and  Dimethyl  Phosphorofluoridate* 


r- 

ON 

vo 

vo 

«n 

tn 

o 

m 

>n 

O 

r- 

r-H 

00 

00 

r-H 

°r 

< 

vo 

00 

On 

ON 

<N 

in 

ON 

CN 

o 

ON 

Tj- 

Tt 

vo 

tn 

p 

T— < 

tri 

r—H 

T— 1 

© 

K 

00 

00 

rf 

r- 

00 

ON 

o 

o 

CN 

CN 

o 

m 

CN 

m 

t> 

r- 

t> 

t> 

t 

VO 

tn 

& 

o 

r-H 

*— 1 

r-H 

r-H 

r-H 

r-H 

^H 

•—< 

•— ' 

r-H 

r-H 

1 

r— H 

i 

r-H 

"N 

rl* 

Uh 

CN 

o 

m 
i n 

o 

o 

cq 

tr- 

oq 

r— H 

O'- 

O 

o 

00 

CN 

00 

r- 

o 

00 

rr 

r-H 

Tf 

ON 

00 

CN 

VO 

tn 

tn 

r-H 

CN 

vo 

r-H 

VO 

Tj- 

< 

cd 

rd 

r-H 

cm’ 

o 

cn 

o 6 

00 

vd 

vd 

t> 

tn 

vd 

vd 

ON 

ON 

O 

r-H 

•ON 

T*H 

r-H 

r-H 

r-H 

ON 

ON 

ON 

00 

00 

00 

00 

CO 

00 

vo 

vo 

c- 

fO 

m 

o 

r— H 

Gh 

m 

o 

cn 

ON 

On 

CN 

On 

VO 

cn 

IP 

0" 

O 

m 

«n 

o 

in 

m 

vo 

m 

On 

tn 

rq 

vo 

r- 

p 

VO 

fq 

vo 

n 

m 

00 

CN 

CN 

$ 

o 

tn 

£ 

00 

to 

t> 

fO 

00 

rd 

CN 

CN 

ON 

^H 

On 

r-H 

On 

On 

00 

^H 

On 

CN 

ON 

o 

o 

CN 

O 

r^ 

r— H 

o 

r-H 

r-H 

r— ( 

rj* 

^H 

r-H 

CN 

r— H 

r-H 

rd 

r— H 

r-H 

r-H 

r-H 

r-H 

r-H 

r-H 

r— H 

rf 

o 

o 

o 

•n 

o 

o 

o 

o 

O 

tn 

o 

O 

00 

O 

o 

o 

o 

o 

Ov 

o 

o 

tit 

o 

On 

00 

ON 

00 

vo 

CN 

ON 

c- 

On 

in 

r-H 

VO 

•n 

fN 

00 

On 

c- 

00 

m 

s 

< 

vo 

tn 

tn 

»n 

m 

tn 

VO 

to 

VO 

VO 

vo 

t> 

r- 

CN 

CN 

CN 

tn 

Ph 

1— 1 

T— 1 

r—H 

r-H 

r-H 

r-H 

r-H 

r-H 

*— H 

CN 

CN 

CN 

rn 

rn 

cd 

o 

c 

m 

o 

o 

o 

o 

o 

o 

o 

00 

o 

o 

r-H 

o 

O 

o 

o 

o 

r-H 

o 

o 

cn 

vo 

CN 

o 

ON 

fN 

i— H 

00 

in 

m 

m 

o 

ON 

CN 

ON 

ON 

m 

CN 

r-H 

m 

00 

00 

»“* 

CN 

» 

l> 

vo 

vo 

vo 

vo 

VO 

tn 

vo 

tn 

*n 

vo 

VO 

VO 

a, 

cd 

rd 

rd 

CN 

cn 

CN 

r-H 

r— H 

r-H 

r-H 

r-H 

r-H 

r— H 

r-H 

r-H 

H 

p 

to 

«> 

r- 

tn 

o 

On 

r- 

tn 

(N 

3 

m 

ON 

m 

s 

VO 

r-H 

tn 

O 

o 

VO 

VO 

0- 

vo 

tn 

to 

Tf 

<N 

(N 

CN 

fN 

r—l 

r- 

VO 

ON 

00 

00 

00 

c- 

TT 

rt- 

>n 

to 

to 

vo 

vo 

vo 

vo 

vo 

vo 

tn 

*n 

tn 

■*fr 

a, 

r-H 

r-H 

i-H 

r— H 

r-H 

r-H 

r-H 

r-H 

r-H 

r-H 

/-“*S 

O 

o 

o 

o 

o 

o 

o 

VO 

O 

CN 

ON 

m 

^r 

ITi 

CN 

m 

r-H 

On 

rn 

vo 

fN 

00 

rn 

CO 

oq 

CN 

»n 

in 

CN 

m 

On 

2.46 

r- 

vo 

tn 

c- 

tn 

uf 

o 

o 

o 

CN 

CN 

cd 

rd 

m 

K 

m 

CN 

tt 

(N 

o 

fN 

O 

fN 

vd 

r-H 

ON 

CN 

CN 

CN 

vd 

CN 

O 

m 

,"H 

CN 

tn 

tn 

tn 

/*“N 

VO 

vo 

o 

vo 

•'fr 

m 

O 

O 

o 

vo 

tn 

rq 

o 

tn 

>n 

o 

ON 

fN 

tn 

o 

»n 

fN 

t- 

p 

o 

ON 

m 

o 

o 

tn 

vo 

O 

O 

r-H 

O) 

00 

oo 

CN 

t> 

w 

tn 

vo 

o 

o 

r-H 

On 

fN 

cn 

ki 

K 

f o 

t-' 

tn 

«n 

*— H 

CN 

Tj* 

ON 

CN 

tn 

CN 

o 

o 

o 

o' 

On 

o’ 

r-H 

m 

-'t 

00 

m 

vo 

VO 

«n 

O 

00 

tn 

o 

w 

m 

00 

tn 

o 

o 

o 

00 

rq 

(N 

oq 

ON 

On 

VO 

oq 

p 

T— H 

ob 

o 

O) 

oq 

°9 

VO 

CO 

o 

o 

r-H 

m 

rd 

m 

00 

m 

tn 

rn 

rn 

O 

CN 

tn 

ON 

CN 

CN 

»n 

CN 

o 

o 

o 

CN 

r-H 

r-H 

CN 

n 

r 

05 

05 

r-i 

CN 

H 

tO 

H 

00 

On 

22 


*The  values  shown  for  each  parameter  are  optimized  using  the  B3LYP/6-3 1 1+G(2d,2p)  level  of  theory  for  the 
phosphinate,  phosphonate,  and  phosphate,  respectively.  The  relative  energies  (in  kcal/mol)  are  relative  to  the 
lowest  energy  structure.  E(z)  and  E(G)  are  the  zero-point  and  Gibbs  free  energies,  respectively.  Bond  lengths 
are  in  angstroms  and  angles  in  degrees. 


4.3 


Neutral  Reaction  Pathway  -  6-Membered  Transition  State. 


To  extend  the  investigation  of  the  hydrolysis  of  the  three  OP  compounds  via 
explicit  water  molecule  solvation,  a  new  model  using  two  water  molecules  was  employed.  This 
methodology  uses  a  less  strained  6-membered  TS  rather  than  the  prior  4-membered  TS.  A 
schematic  for  the  two-water  neutral  reaction  pathway  is  depicted  in  Scheme  3  with  the 
significant  atoms  numbered.  The  nucleophile  03  oxygen  approaches  the  phosphorus  atom  anti 
to  the  fluorine-leaving  group  as  in  the  prior  anionic  and  one-water  reaction  pathway.  However, 
this  mechanism  transfers  a  proton  to  a  second  water  molecule,  which  in  turn  concurrently 
transfers  its  proton  to  the  phosphoryl  oxygen  to  stabilize  the  negative  charge.  This  produces  a 
unique  6-membered  TS  with  dual  proton  transfer  occurring  simultaneously  with  the  PI -03  bond 
formation,  resulting  in  six  bonds  being  made/broken  concurrently  (counting  the  double  to  single 
phosphoryl  bond)  as  is  evidenced  by  the  single  imaginary  frequency.  The  H5  proton  then 
departs  via  a  similar  6-membered  TS  with  a  water  and  the  fluoride-leaving  group  forming  the 
phosphate,  water  and  hydrogen  fluoride. 

The  structural  and  energetic  results  for  the  complete  gas  phase  two-water  neutral 
reaction  pathways  for  the  hydrolysis  of  dimethylphosphinic  fluoride,  methyl 
methylphosphonofluoridate,  and  dimethyl  phosphorofluoridate  are  consolidated  in  Table  6.  For 
simplicity  in  discussing  results  in  the  text  of  the  two-water  pathway,  the  first,  second  and  third 
results  reported  are  of  the  phosphinate,  phosphonate  and  phosphate,  respectively.  In  Figure  7, 
the  PES  of  the  phosphonate  is  displayed,  together  with  the  corresponding  gas-phase  optimized 
structures  depicted  relative  to  their  positions  along  the  reaction  coordinate.  The  PESs  of  the 
phosphinate  and  phosphate  are  similar  to  the  phosphonate  PES,  and  thus,  are  not  depicted 
separately.  The  reactants  and  products  at  infinite  separation  are  not  depicted. 

As  with  the  one-water  hydrolysis,  only  the  lowest  energy  conformations  are 
reported.  With  the  two-water  pathway,  there  are  several  possible  conformations.  For  example, 
as  with  the  single  water  analog,  there  are  viable  conformations  with  the  methoxy  group(s)  on  the 
same  side  as  the  water(s)  as  well  as  on  the  side  with  the  fluoride-leaving  group.  With  each  of 
these  conformations,  the  water  protons  not  in  the  TS  ring  give  rise  to  two  and  four  distinct 
conformations  for  the  phosphonate  and  phosphate,  respectively. 

It  is  again  noticed,  as  in  the  one-water  hydrolysis,  that  the  phosphate  has 
tighter/stronger  bonds  throughout  the  PES  thereby  pointing  to  the  significance  of  substituents  on 
the  reaction  rates.  Also,  TB  structure  4  is  more  stable  than  TB  structure  6  only  for  the  phosphate 
by  2  kcal/mol,  which  may  contribute  to  the  phosphate’s  lower  reactivity.'  A  major  difference 
between  the  one- water  and  two-water  systems  was  in  not  locating  HB  structure  8  for  the  two- 
water  phosphate  system.  By  comparing  the  three  OP  systems  at  the  B3LYP/6-3 1 1+G(2d,2p) 
level  of  theory,  the  free  energy  barriers  for  the  two-water  hydrolysis  are  32.1,  34.4  and  37.0 
kcal/mol,  respectively.  This  correlates  to  a  small  1.5,  2.0  and  5.1  kcal/mol  reduction  over  the 
strained  one-water  hydrolysis  values.  Therefore,  it  is  feasible  that  the  one-  and  two-water 
mechanisms  are  possible. 
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Table  6.  Structural  and  Energetic  Results  for  the  Two- Water  Hydrolysis  of  Dimethylphosphinic  Fluoride, 
Methyl  Methylphosphonofluoridate,  and  Dimethyl  Phosphorofluoridate* 
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*The  values  shown  for  each  parameter  are  optimized  using  the  B3LYP/6-3 1 1+G(2d,2p)  level  of  theory  for  the 
phosphinate,  phosphonate,  and  phosphate,  respectively.  The  relative  energies  (in  kcal/mol)  are  relative  to  the 
lowest  energy  structure.  E(z)  and  E(G)  are  the  zero-point  and  Gibbs  free  energies,  respectively.  Bond  lengths 
are  in  angstroms  and  angles  in  degrees. 


kcal/mot 


Figure  7.  Relative  Free  Energies  for  the  Two-Water  Neutral  Reaction  Pathway 
for  the  Hydrolysis  of  Methyl  Methylphosphonofluoridate 


Plot  of  the  relative  free  energies  for  the  two-water  neutral  reaction  pathway  for  the  hydrolysis  of  methyl 
methylphosphonofluoridate.  Gas  phase  optimized  geometries  are  shown  relative  to  their  position  along  the 
reaction  coordinate.  Methyl  hydrogens  omitted  for  clarity. 


4.4  Pathways  Syn  to  the  Fluorine  Atom. 

It  is  possible,  though  unlikely,  that  the  hydrolysis  path  involves  nucleophilic 
attack  on  the  fluorine  side  of  the  molecule  rather  than  backside  attack  and  thereby  eliminates  the 
TB  structure. 


The  corresponding  one-  and  two-water  TSs  were  determined  using  the  B3LYP/6- 
3 1 1  +G(2d,2p)  level  of  theory  with  the  phosphinate.  The  one-water  barrier  for  the  phosphinate  is 
20.02  kcal/mol  higher  than  the  corresponding  barrier  with  the  water  anti  to  the  fluorine-leaving 
group.  The  two-water  barrier  for  the  phosphinate  is  25.53  kcal/mol  higher  than  the 
corresponding  barrier  with  the  water  anti  to  the  fluorine-leaving  group.  Therefore,  these 
pathways  do  not  seem  feasible. 
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5. 


CONCLUSIONS 


The  hydrolysis  of  three  OP  compounds  —  a  phosphinate,  a  phosphonate,  and  a 
phosphate  —  was  investigated  in  the  gas  phase  and  in  solution  using  a  hydroxide  anion  (anionic 
reaction  pathway)  and  a  water  molecule  (neutral  reaction  pathway).  The  instability  of  hydroxide 
anions  in  the  gas  phase  led  to  transition  states  that  were  lower  in  energy  than  both  the  reactants 
and  products;  therefore,  solvation  was  needed  to  accurately  describe  the  anionic  reaction 
pathway.  The  SCI-PCM  SCRF  solvation  model  implemented  in  Gaussian  94  was  used  to 
investigate  the  anionic  hydrolysis.  The  electronic  energy  profile  was  compared  to  the  free 
energy  profile  for  the  gas  phase  and  SCI-PCM  methods.  The  free  energy  profile  gives 
insignificant  differences  when  compared  to  the  electronic  energy  profile  for  the  gas  phase 
calculations  and  only  needs  to  be  taken  into  account  with  the  SCI-PCM  calculations.  The  SCI- 
PCM  solvation  model  produced  an  initial  energy  barrier  with  the  incoming  nucleophile  and  had 
an  overall  downhill  path  to  produce  products  for  all  three  systems.  The  existence  probability  of 
the  pentacoordinate  species  (which  is  a  stable  intermediate  in  gas  phase  calculations)  is  greatly  - 
reduced  with  the  inclusion  of  solvation.  The  correct  qualitative  trends  in  the  reaction  rates  of  the 
three  OP  compounds  were  produced  with  the  solvation  method;  however,  other  solvation 
methods  are  being  evaluated  for  greater  quantitative  precision. 

Due  to  the  high  cost  associated  with  the  SCI-PCM  solvation  calculations  (i.e., 
numerical  second  derivatives),  a  neutral  reaction  pathway  methodology  was  employed  to 
generate  quantitative  results  at  a  reasonable  cost.  This  methodology  is  based  upon  using  a  single 
proton  from  the  nucleophile  to  stabilize  the  charge  on  the  system.  The  complete  neutral  reaction 
pathway  of  the  phosphinate  at  the  RHF/6-31  lG(2d,2p),  B3LYP/6-31 1+G(2d,2p),  and  MP2/6- 
31 1+G(2d,2p)  levels  of  theories  showed  an  initial  rate  determining  step  with  the  attack  of  the 
incoming  nucleophile.  The  B3LYP  hybrid  functional  gives  structural  and  energetic  results  of 
similar  accuracy  as  the  more  expensive  MP2  level  of  theory  using  the  same  basis.  This  neutral 
reaction  pathway  methodology  using  the  B3LYP  level  of  theory  is  currently  being  used  to 
investigate  the  phosphorylation  of  the  active  site  of  acetylcholinesterase. 

Solvation  calculations  using  30(+)  water  molecules  are  currently  being 
investigated  using  the  B3LYP  level  of  theory  (to  be  reported  elsewhere).  As  computational 
power  has  become  readily  available,  these  calculations  have  become  viable  and  are  very  similar 
in  cost  to  the  expensive  SCRF  SCI-PCM  solvation  method,  without  using  the  empirical 
assumptions  associated  with  such  methods. 

6.  APPLICATIONS 

As  mentioned  earlier,  computational  chemistry  is  a  tool  for  supporting  analytical 
and  developmental  efforts.  This  study  characterizing  the  transition  structures  and  the 
corresponding  energies  for  the  hydrolysis  of  phosphorus  esters  has  application  to  the  following 
topics: 


a.  Analysis  and  evaluation  of  interaction  of  novel  chemical  threat  agents 
with  acetylcholine  esterase  and  other  biological  receptors.  The  principal  value  will  be 
determination  of  the  energies  of  activation  for  inhibition  reactions  to  determine  whether  the 
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novel  agents  are  sufficiently  reactive  (naturally,  as  well  as  catalyzed)  to  inhibit  the  biological 
process. 


b.  Determination  of  the  energies  of  activation  for  reactions  of  CW  agents 
with  nucleophiles,  oxidants,  and  other  potential  decontaminants.  The  chemical  reactions  of 
candidate  decontaminants  with  traditional  chemical  warfare  agents  in  a  variety  of  solvents  will 
be  determined  to  predict  the  products  and  the  reaction  rates  at  militarily  significant  temperatures. 
These  calculations  will  help  identity  the  candidates  whose  reactivity  falls  within  the  desired 
range  that  is  sufficiently  reactive  to  destroy  the  agents  but  mild  enough  to  minimize  hazard  to 
personnel,  equipment,  and  the  environment. 

c.  Persistency  of  agents  in  the  environment.  The  length  of  time  that  an  agent 
remains  in  the  environment  depends  upon  its  vapor  pressure  and  resistance  to  hydrolysis  in 
conjunction  with  environmental  conditions  such  as  temperature  and  precipitation,  etc.  Results 
from  these  calculations  on  the  individual  molecule  provide  a  starting  point  for  future  work. 
Creating  a  mathematical  surface  that  represents  soil  particles  and  determining  the  effects  of 
molecular  interactions  is  not  a  trivial  exercise;  however,  earlier  semiempirical  calculations 

24.25 

indicated  the  approach  was  feasible. 
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